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GMSK 十 PN 信号 信 噪 比 估计 
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摘 要 : 针对 数 传 、 测 距 一 体 化 的 调制 方式 高 斯 最 小 频 移 键 控 复合 伪 码 测 距 (Gaussian minimum shift keying 十 
pseudo-noise，GMSK 十 PN) 信 号 的 信 噪 比 估计 ,提出 了 一 种 精度 较 高 复杂 度 较 低 的 快速 算法 。 该 算法 采用 层 状 
结构 分 解 的 方式 ,将 信号 分 解 为 若干 信号 单元 ,对 每 一 个 分 解 的 信号 单元 进行 功率 分 配 计算 或 信 噪 比 估 计 , 最 后 
快速 准确 地 计算 出 整个 信号 的 信 噪 比 。 信 噪 比 在 1~13 dB 之 间 时 ,估计 误差 小 于 0.5 dB。 利 用 Simulink 仿真 工 
具 , 验 证 了 该 算法 的 有 效 性 。 
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SNR estimation algorithm for GMSK 十 PN signal 


SHI Liguo'’, WANG Zhugang'’, XIONG Weirmingl 
(1. Center for Space Science and Applied Research ,Chinese Academy of Sciences, Beijing 100190, China; 
2. University of Chinese Academy of Sciences, Beijing 100190, China) 


Abstract: Gaussian minimum shift keying 十 pseudornoise (GMSK 十 PN) is a modulation method which a- 
chieves Simultaneously transmission of GMSK data and PN. For the signal to noise ratio (SNR) estimation of 
GMSK 十 PN modulation, a fast algorithm with high accuracy and low complexity is proposed. This algorithm 
applies layered structure decomposition to decompose signals into a number of units. Each decomposed signal’s 
power is calculated or SNR is estimated. Then, signal power distribution calculation or SNR estimation is con- 
ducted successively for according units. Finally, SNR of the entire signal is calculated quickly and accurately. 


Estimation error is less than 0.5 dB when SNR is between 1~13 dB. Simulink simulation results verify the ef- 
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fectiveness of the proposed algorithm. 
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分 布 式 卫星 系统 要 实现 编队 飞行 ,需要 在 测定 子 卫 星 
之 间 相 对 位 置 的 同时 ,还 要 求 有 子 卫星 间 的 数据 交互 。 测 
距 \ 通 信 一 体 化 设计 在 许多 卫星 通信 系统 中 都 有 涉及 ,例如 
中 国 科学 院 空间 科学 与 应 用 中 心 的 深 空 应 答 机 。 但 其 采用 
的 一 体 化 设计 方案 是 时 分 复 用 的 方式 ,无 法 实现 同时 完成 
测 距 和 通信 两 种 任务 。 目 前 ,空间 数据 系统 咨询 委员 会 
(consultative committee for space data systems,CCSDS) 正 
在 酝酿 一 种 可 以 同时 传输 高 码 率 数 传 信号 和 测 距 信号 的 技 
术 体 制 。 它 采用 高 斯 最 小 频 移 键 控 (Gaussian minimum 
shift keying,GMSK) 进行 数 传 信号 的 传输 ,采用 再 生 伪 


ll 


码 (pseudo-noise,，PN) 进 行 测 距 ” ,被 称 为 高 斯 最 小 频 移 键 
控 复 合 伪 码 测 距 (Gaussian minimum shift keying 十 pseudo- 
noise,GMSK 十 PN) 调 制 方式 ” ,并 成 为 NASA 月 球 任务 的 
备 选 方案 之 一 。 文 献 L8] 最 早 对 GMSK 十 PN 调制 方式 
的 测 距 信号 副 载波 实现 形式 进行 了 频谱 分 析 , 并 给 出 了 仿 
真 结 果 。 文 献 L9] 研 究 了 更 具 一 般 性 的 连续 相位 调制 (con- 
tinuous phase modulation,CPM) 下 的 测 距 一 体 化 设计 调制 
解 调 算法 。 文 献 [L10-11] 给 出 了 GMSK 十 PN 和 非 平衡 四 相 
相 移 键 控 (unbalanced quaternary phase shift keying， 
UQPSK) 的 优 缺 点 和 性 能 差异 。 文 献 [12-13 | 对 GMSK 十 
PN 无 副 载波 实现 形式 的 相位 抖动 和 同步 进行 了 仿真 。 
GMSK 十 PN 调制 方式 从 信号 层面 上 来 看 是 真正 的 测 距 通 
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信 一 体 化 , 它 兼 有 GMSK 和 PN 测 距 的 优点 29 : 旁 办 衰减 
快 、. 频 谱 紧 凑 、 带 宽 效率 高 .适合 饱和 功放 ; 测 距 捕获 时 间 
短 ,硬件 实现 复杂 度 低 , 测 距 精度 相对 较 高 。 由 于 采用 同一 
个 链 路 同时 完成 数 传 和 测 距 两 个 功能 ,这 种 调制 方式 也 有 
相应 的 缺点 : 解 调 算法 相对 复杂 , 数 传 信号 和 测 距 信 号 互 为 
噪声 ,影响 了 解 调 的 性 能 。 因 此 ,为 了 保证 数 传 信号 的 误 码 
率 性 能 和 测 距 的 精度 , 数 传 信号 采用 相干 解 调 的 方式 。 信 
噪 比 估计 是 信道 估计 的 一 个 重要 组 成 部 分 , 它 提供 了 切换 、 
功率 控制 和 信道 分 配 算 法 所 需 的 信道 质量 信息 。 在 解 调 过 
程 中 ,许多 解 调 和 译 码 的 算法 都 需要 准确 的 信 噪 比 信息 才能 
达到 最 优 的 性 能 "…] 。 本 文 根 据 GMSK 十 PN 信号 的 特点 ， 
提出 了 一 种 适用 于 GMSK 十 PN 信号 的 信 噪 比 估 计算 法 。 


1 GMSK 十 PN 调制 解 调 技术 


1.1 调制 技术 
基于 GMSK 十 PN 的 调制 体制 能够 同时 进行 数 传 和 
测 距 的 系统 框图 如 图 1 所 示 。 


图 1 GMSK 十 PN 信号 调制 框图 


测 距 符号 经 过 电 平 转换 和 半 正 弦 滤 波 后 , 乘 以 PN 测 
距 调 制度 (也 称 测 距 调 制度 ) 得 到 PN 测 距 信号 相位 ; 数 传 
符号 经 过 预 编码 、 电 平 转换 、 高 斯 滤波 和 积分 器 后 ,得 到 
GMSK 相位 ;PN 相位 和 GMSK 相位 相 加 后 再 经 过 相位 调 
制 , 即 可 得 到 GMSK 十 PN 调制 信号 
s(t) = V2P cos(wt+ Bowsk (DT mapes(t)) (1) 
其 中 


prs (D) = DD sin( 3 ort) (2) 


式 中 ,Pu 为 信号 的 总 功率 ;w. 为 载波 角 频 率 ; pevsk (7?) 为 数 传 
信号 调制 相位 ;pm () 为 PN 测 距 信号 相位 ; DC) 为 PN 序列; 
m4 为 PN 测 距 调制 度 ;wm 为 PN 序列 角 频 率 。 
1.2 解 调 技术 

因为 基带 的 GMSK 数 传 信号 (所 代表 的 相位 ) 和 PN 测 
距 信号 (所 代表 的 相位 ) 在 时 域 上 是 按 比例 相 加 的 ,所 以 如 
果 得 到 其 中 的 任何 一 种 信号 ,就 可 以 由 复合 信号 与 某 一 种 
“再 生 ” 的 信号 “ 相 减 ”, 得 到 另 一 种 信号 。 因 为 按照 常规 ， 
GMSK 信号 在 复合 信号 中 所 占 的 能 量 比例 较 大 。 所 以 ,可 
以 直接 在 已 调 载波 信号 上 解 调 ,得 到 GMSK 基带 信号 。 再 
由 GMSK 基带 信号 的 再 调制 得 到 的 相位 ,与 输入 信号 进行 
复 相 关 , 即 可 得 到 单独 的 PN 测 距 信号 525 ( 见 图 2) 。 


言 号 悄 


GMSK 再 调制 GMSK 信 号 剥离 


图 2 GMSK 十 PN 信号 解 调 框图 
2 GMSK 十 PN 信号 信 噪 比 估 计 


2.1 信号 简化 波形 

由 于 GMSK 十 PN 信号 既 包 括 数 传 信号 ,又 包括 伪 码 
测 距 信号 ,无 法 直接 进行 信 噪 比 估 计 。 而 GMSK 十 PN 信 
号 的 结构 非常 特殊 ,在 数 传 信 号 的 基础 上 ,直接 加 入 了 PN 
测 距 信号 (作为 数 传 信 号 的 相位 抖动 存在 ) , 故 可 以 采用 逐 
层 简化 的 方式 ,将 GMSK 十 PN 信号 分 解 为 若干 层 , 如 图 3 
所 示 。 本 文采 取 反 推 的 方法 , 先 估计 出 PN 信号 中 C1 子 码 
分 量 的 信 噪 比 , 即 可 以 通过 PN 信号 与 C1 子 码 分 量 的 功率 
分 配 ,得 到 PN 信号 的 信 噪 比 。 进 而 ,通过 GMSK 十 PN 信 
号 的 能 量 分 配 ,得 到 整个 GMSK 十 PN 信号 的 信 噪 比 。 


Ma 


KN 201g() 


图 3 信 噪 比 估计 简化 模型 


在 图 3 中 ,第 1 层 ( 最 内 圈 ) 的 面积 代表 噪声 的 功率 
Nm ,第 2 层 ( 最 内 轿 和 次 内 圈 之 间 ) 的 面积 代表 PN 测 距 信 
号 C1 子 码 分 量 的 功率 Pu ,第 3 层 (次 内 圈 和 次 外 圈 之 间 ) 
的 面积 代表 PN 测 距 信号 C2~C6 子 码 分 量 的 功率 Pos_c， 
第 4 层 ( 次 外 层 和 最 外 层 之 间 ) 的 面积 代表 数 传 GMSK 信 
号 的 功率 Pewsk 。Pa /Nu 可 以 通过 最 常用 的 均值 方差 比 可 
以 求 得 ; 当 测 距 PN 的 码 型 确定 之 后 , Pex/Poa 是 确定 的 ; 
Pousk /Pex 可 以 通过 mv 来 计算 。 简 化 模型 如 图 3 所 示 。 这 
种 新 的 信 噪 比 估计 算法 ,可 以 大 大 降低 GMSK 十 PN 信号 
的 信 噪 比 计 算 的 复杂 度 。 
2.2 信 噪 比 估计 
2.2.1 数 传 . 测 距 信号 能 量 分 配 

由 贝 塞 尔 函 数 可 知 

sin[ msin(at ) ] 一 


2>) Ti Cn) sinCkat) 2 Cm) sin(at) (3) 
kodd 
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cos[Lmsin(at ) | 中 Cl 分 量 成 型 后 即 为 正弦 波 ， PN 信号 成 型 后 波形 如 图 5 
Jo On) 42D Cm eoshat) SS Jom) (4) ”所 示 。 
式 中 ,J (m) 为 第 一 类 k& 阶 贝 塞 尔 画 数 。 
GMSK 十 PN 的 调制 信号 表达 式 如 式 (1) 所 示 , 将 式 (3)、 ge 
式 (4) 代 入 式 (1) 并 将 其 展开 车 
s(1) = V2Pio Jo ma)sinCwt 十 peusk (7)) 四 显 下 
V2Pit msin( (3 二 wpn 十 ww )* +t gousx(D ) 十 ee 
部 (a) 半 正 弦 成 型 Tausworhte 码 
2Pm 1 Cma) sin( (Zon We 有 一 Yovsk(D】 (5) | 
对 于 数 传 信 号 的 解 调 ,所 用 的 是 式 (5) 中 等 号 右边 的 第 
一 项 ;对 于 测 距 信号 的 解 调 , 所 用 的 是 式 (5) 中 等 号 右边 的 
第 2 项 和 第 3 项 。 因 此 , 今 数 传 信 号 功率 为 Pewsk ,PN 测 距 
信号 的 功率 为 Pa , 则 有 t/ms 
(b) 半 正 弦 成 型 Tausworhte 码 C1 子 码 分 量 
Pems 2 
PB = Jo (ma) 图 5 ” Tausworthe 码 成 型 波形 
Ei (6) 经 过 整形 后 的 PN 与 正弦 波 相 比 ,其 功率 有 部 分 衰减 
i 册 式 (8) 记 示 。 
PIN 23 (ma) PN 
Poi a 10lg( BE) 一 一 20lg(6 ) (8) 


式 中 ,ww 为 测 距 调制 度 , 是 GMSK 十 PN 信号 中 GMSK 信 
号 和 PN 信号 能 量 分 配 的 唯一 参数 。 随 mm 变化 ,GMSK 十 
PN 信号 及 其 数 传 , 测 距 部 分 能 量 的 归 一 化 比例 的 变化 趋势 
如 图 4 所 示 。 


0.8 oe hi tt 
0. 
祥 0.6 bbs ee, Sh isa pip eis bt Wyte “ie 
访 05 
不 0. 针 --1--- 
i a 


0 0.1 505 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
测 距 调 制度 /rad 
-6-: GMSK 功 率 分 配 ; -6 : PN 信号 功率 分 配 ; 
-6-: 总 功率 衰减 误差 ;4-: GMSK/PN 功 率 分 配 。 


图 4 GMSK 与 PN 信号 功率 分 配 


当 wm 为 确定 值 时 ,Povsr /Prx 也 是 一 个 固定 值 。wmw 的 
值 偏 低 时 (小 于 0. 2 rad) , 则 测 距 部 分 能 量 较 小 ,影响 测 距 
情 度 ;wm 的 值 偏 高 时 (大 于 0. 45 rad) , 则 数 传 部 分 能 量 较 
小 ,影响 数 传 解 调 性 能 "”*” 。 因 此 ,要 合理 选择 ws 。 

由 于 Bessel 画 数 存在 近似 , 随 着 ma 的 变化 ,GMSK 信 
号 的 总 功率 P, 存 在 不 同 程度 的 衰减 。 所 以 在 求 得 数 传 信 
号 功率 Povsrk 和 PN 测 距 信号 功率 Pv 后 ,需要 对 GMSK 十 
PN 信号 总 功率 Pi 进行 修正 ,如 式 (7) 所 示 。 


PP 一 Pewsk 十 Ppx 万 Ca ) PN (7) 
六 Gn) +2 1m) 2 GM) JGnm)+l1 
2.2.2 PN 测 距 信号 功率 分 配 


PN 测 距 信号 经 过 半 正 弦 成 型 后 ,与 正弦 波 极 其 相似 ,其 


式 中 ,5 为 Tausworthe 码 中 C1 分 量 与 复合 码 同 相 的 比例 ， 
其 值 随 码 型 变化 。 可 以 求 得 
Pa _ 1 
本 一 宇 
对 于 平衡 2 倍加 权 Tausworthe 码 (weighted-2 balanced 
Tausworthe ,T2B) ,有 


(9) 


Ppyn 
Po 

对 于 平衡 4 倍加 权 tausworthe 码 (weighted-4 balanced 
Tausworthe ,T4B) ,有 


10lg gf ) 一 4 049 6 dB 


10lg( PE 
2.2.3 单 音信 号 信 噪 比 估计 

PN 测 距 信和 号 的 C1 分 量 ,经 过 半 正 弦 成 型 滤波 后 , 即 为 
正弦 波 ,而 对 这 种 单 载波 恒 包 络 的 信号 的 信 噪 比 估 计 可 以 
使 用 最 简单 的 均值 方差 比 ”" 求 得 
SS. 101 | ElaCG) 了 | 


20549 7 .dB 


N VarlL a(n) | (10) 


式 中 ,a() 为 PN 测 距 信号 C1 分量 瞬时 包 络 。C1 分 量 采 
样 模型 如 图 6 所 示 。 其 中 ,i, 为 采样 间隔 ;i 为 测 距 码 片 周 
期 ;A 为 实际 接收 信号 的 幅度 


ts 
1 


图 6 PN 半 正 弦 成 型 信号 采样 量化 模型 
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使 A' 为 信号 幅度 的 估计 值 ,x(k) 为 每 个 采样 点 的 幅 
度 , 则 ,经 过 采样 后 估计 的 chip 幅度 为 


Dr) 天 | zcod i Asin( TL) dt 


A’” 2 0 a 
t, 40 t, ts 
于 Atcos( 开 )| = 对 (11) 
A te 0 nt, 
信号 功率 的 估计 值 S 为 
i ntA’” 外 A 省 2 / 
S 2 。( 21. ) 5 (22 ) EFE(A) (12) 
噪声 功率 估计 值 N 为 
N = EVar(A’) (13) 


式 中 ,C1 分 量 理 论 信 噪 比值 为 总 -, 其 信 噪 比 估计 值 司 


为 


SM 


--------+--------+--- 


ee SR 


4 6 8 10 12 14 
信 品 比 实际 值 /dB 

一 一 : 实际 值 ; -所 -: 估计 值 。 

(a) 信 品 比 实际 值 与 估计 值 关 系 


经 过 仿真 可 以 发 现 , 信 噪 比 估计 值 与 实际 值 非常 吻合 ， 
误差 在 0.5 dB 以 下 。 参 与 计算 的 码 片 数 在 大 于 1 000 时 ， 
信 噪 比 估计 值 基本 稳定 实际 值 附近 。 


4 结 论 


本 文 给 出 了 在 加 性 高 斯 白 噪声 信道 中 基带 GMSK 十 
PN 信号 的 信号 模型 ,同时 给 出 了 在 此 模型 下 的 简单 ,快捷 
的 信 噪 比 算法 。 并 通过 Simulink 对 该 算法 进行 了 仿真 。 该 
算法 精确 度 较 高 ,计算 复杂 度 低 , 且 参 与 计算 的 码 片 数 目 较 
少 , 具 有 工程 的 可 实现 性 。 
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图 7 仿 


Pa Bs ECA) 
N, NN’, 8 1t. Var(A’) 2f. 
由 第 2.2.1 节 、 第 2.2.2 节 和 第 2.2.3 节 的 结论 ,可 以 


得 到 GMSK 十 PN 信号 总 功率 与 噪声 功率 之 比 为 


(14) 


Pe J Ona) . 

No 2J2 (ma)Ji (ma)+t+l 

1 1 xt EC(A)  A = 
六 人 


3 结果 分 析 


通过 Simulink 搭建 仿真 模型 , 取 采 样 率 为 f. 二 100 MHz， 
PN 速率 为 f. 二 2 MHz, 测 距 调 制度 wj 二 0.2 rad, 再 生 PN 
采用 T4B, 测 得 信 噪 比 估计 值 如 图 7 所 示 。 
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